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Numerical models of rotating neutron stars are constructed for four equations of state using the 
computer code rns written by Stergioulas. For five selected values of the star's gravitational mass 
(in the interval between 1.0 and 1.8 solar masses) and for each equation of state, the star's angular 
momentum is varied from J = to the Keplerian limit J = Jmax- For each neutron-star configuration 
we compute Q, the quadrupole moment of the mass distribution. We show that for given values of 
M and J, \Q\ increases with the stifi'ness of the equation of state. For fixed mass and equation of 
state, the dependence on J is well reproduced with a simple quadratic fit, Q ~ —aJ^jMt?^ where c 
is the speed of light, and a is a parameter of order unity depending on the mass and the equation 
of state. 

Subject headings: stars: neutron — stars: rotation 
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I. INTRODUCTION 

A rapidly rotating neutron star is not spherically sym- 
metric, because the rotation creates a deformation in the 
stellar mass distribution. In turn, the star's oblateness 
creates a distortion in the gravitational field outside the 
star, which is measured by Qab, the quadrupole- moment 
tensor. Our purpose in this paper is to compute the 
quadrupole moments of rotating neutron stars, under the 
assumptions of rigid rotation, axial symmetry, and reflec- 
tion symmetry about the equatorial plane. 

In Neviftonian theory, the gravitational field outside 
any nonspherical body is given by (Jackson 1975) 
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where M is the body's mass and r = \x\ = {SabXaXbY^'^ 
is the distance from the center of mass; summation over 
repeated indices is understood. In terms of e, the mass 
density inside the body, the quadrupole-moment tensor 
is given by 



ab= e{x') {x'^x'f, - Uabr''^)d^x'. 



(2) 



Thus, in Newtonian theory it is possible to define the 
quadrupole-moment tensor either in terms of the falloff 
behavior of the gravitational potential, or in terms of an 
integral over the mass density. 

If the body is axially symmetric about the z direction, 
then Eqs. (0) and @) reduce to 



and 



Qab = diag(--i(3, -iQ, f Q), 



(3) 



(4) 



where cos9 = z/r, P2{x) = i(3a;^ — 1), and 

Q= [ e{r',e')r'^ P2{cos9')d^x'. (5) 



If the body is also reflection symmetric about the equa- 
torial plane {9 = tt/2), then the 0{l/r*) term vanishes 
identically and Eq. (^) is valid up to terms of order 1/r^. 

It is convenient to introduce the dimensionless quantity 
q, related to Q and M by 



G2 M3' 



(6) 



If we also introduce a dimensionless potential $* = $/c^ 
and a characteristic length M* = GM/c^, then Eq. (||) 
becomes $*(r,6l) = -M*/r - q{M* /rfP2{cos9) + 
0{l/r^). This equation resembles more closely the rela- 
tivistic analogue displayed in the following section. 

Because neutron stars are compact objects with strong 
internal gravity, their gravitational fields must be de- 
scribed within the framework of general relativity. Equa- 
tions (|l|)-(^ are therefore not valid for neutron stars. In 
Sec. 2 of this paper we describe an operational way, due 
mostly to Fintan Ryan (1995), of defining the quadrupole 
moment of an axially symmetric body in general relativ- 
ity^. In Sec. 3 we describe a method for computing the 
quadrupole moment of a rotating neutron star, under the 
assumptions of rigid rotation, axial symmetry, and re- 
flection symmetry about the equatorial plane. Our com- 
putations rely on a numerical integration of the hydro- 
static and Einstein field equations for selected equations 



*In general relativity, multipole moments are defined for 
both the mass density e and the current density j* — ev, 
where v is the fluid velocity. Throughout this paper, the term 
"quadrupole moment" will refer specifically to the quadrupole 
moment of the mass distribution. 
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of state, carried out with a numerical code written and 
made publicly available by Nikolaos Stergioulas (1995). 
In Sec. 4 we present and discuss our results. 

Computations of quadrupole moments of rotating neu- 
tron stars were presented before in the literature: Datta 
(1988) has carried out such computations for a wide selec- 
tion of equations of states, but his results are restricted 
by an assumption of slow rotation; on the other hand, 
Salgado et al. (1994a and 1994b) have computed, also 
for many equations of state, quadrupole moments for 
maximum-mass configurations of rapidly-rotating neu- 
tron stars. While our own analysis does not incorporate a 
very large selection of equations of state, it explores more 
thoroughly the relationship between quadrupole moment 
and angular momentum, both for slow and rapid rota- 
tions. 

The motivation for this paper comes from the realiza- 
tion that the gravitational waves emitted by a binary 
system of rotating neutron stars will be affected by the 
stars' quadrupole moments. Although this effect is small, 
it is comparable in magnitude to that due to the general 
relativistic spin-spin interaction (Kidder, Will, & Wise- 
man, 1993; Apostolatos et al. 1994). This application 
of our results is discussed in a separate paper (Poisson 
1998). 



II. QUADRUPOLE MOMENTS OF AXIALLY 
SYMMETRIC BODIES IN GENERAL 
RELATIVITY 

Equation (^) states that the quadrupole moment Q is 
the coefficient of the "P2 (cos 6')/r'^" part of the Newto- 
nian potential. As we shall see, a similar statement holds 
in general relativity. However, because of the nonlinear- 
ities of the Einstein field equations, it is more difhcult 
to express the quadrupole moment in terms of an inte- 
gral over the source. We shall come back to this point in 
Sec. 3. 

The metric of a stationary, axially symmetric body can 
be written in the form (Bardeen & Wagoner 1971) 
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where the potentials v, B, uj, and a are functions of r 
and 9. Butterworth and Ipser (1976) have calculated the 
asymptotic behavior of these potentials for large r. Apart 
from a slight change in notation, they find 



V = -M/r + \b(Mlrf 



- q{M/rf P2{cose) + 0{l/r^), 
B = l + b(M/rf + 0{l/r^), 
Muo = 2x{M/rf + 0{l/r^). 



(8) 
(9) 
(10) 



An expression for a will not be needed. Here, M = 
GMjc? is the mass of the body, and x is a dimensionless 
measure of its angular momentum: 



where J is the body's angular momentum; as these equa- 
tions indicate, we work with units such that G — c — 1. 
The factors of G and c will be re-inserted in Sec. 4. 

Equations (||)-(|l^) involve two additional (dimension- 
less) parameters, b and q. The parameter b was left un- 
determined by Butterworth and Ipser (1976). However, 
it is not a free parameter, and it is possible to calculate 
it by taking the spherically symmetric limit of the metric 
(0)-(|o|) and comparing with the Schwarzschild solution 
in isotropic coordinates (Schutz 1985). Since these must 
be identical, we find 



4 



(12) 



On the other hand, g is a free parameter. Because qM^ 
is the coefhcient of the "P2 (cos 6')/r'^" part of v (which 
may be identified as the Newtonian potential) , this is the 
quantity that should be interpreted as the quadrupole 
moment. The parameter g, then, would be identified 
with the dimensionless quantity introduced in Eq. (|^). 

The difficulty with this argument is that it is based 
upon the specific coordinate system of Eq. (^ , and the 
question arises as to how it would be affected by a coordi- 
nate transformation. As we shall see, the statement that 
qM^ is the quadrupole moment of the body is true irre- 
spective of the coordinate system. The rest of this section 
is devoted to a proof of this important statement. 

There exists a well developed literature on coordinate- 
invariant characterizations of multipole moments in gen- 
eral relativity. Such work was pioneered by Geroch 
(1970) and pursued by Hansen (1974) and Thorne (1980). 
Here we follow the simple prescription due to Ryan 
(1995), which is built on previous work by Fodor, Hoense- 
laers, and Perjes (1989). Ryan's prescription relies on 
the existence of a general relationship between the set 
of coordinate-invariant multipole moments defined by 
Hansen (1979), and a coordinate-invariant quantity Ai? 
defined as follows. 

Let a test particle move, in the absence of any exter- 
nal force, around a central body of mass M and angular 
momentum J, in such a way that f2, its angular veloc- 
ity as measured at infinity, is uniform. (Such orbits will 
be termed "circular".) The central body is assumed to 
be stationary, axially symmetric, and reflection symmet- 
ric about the equatorial plane. We assume also that the 
motion of the test particle is confined to the equatorial 
plane, and that the particle moves in the same direction 
as the body's rotation. [If <^ denotes the angle of rota- 
tion about the symmetry axis (which coincides with the 
body's rotation axis) and t is the time measured at in- 
finity, then = d(t)/dt > 0. Such coordinates are used in 
Eq. (0).] 

The orbital energy per unit test-particle mass, denoted 
E, is a coordinate-invariant quantity which depends on 
the orbital motion of the particle. Because each circular 



2 



orbit can be labeled with the value of its angular velocity 
(also a coordinate- invariant quantity), we have that 
E can be expressed as a function of O. The quantity 
considered by Ryan (1995) is 



AE 



(13) 



This evidently is a coordinate-invariant quantity, and 
Ryan derives the following expression for it: 



AE 
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(14) 



where Q (denoted M2 by Ryan) is the coordinate- 
invariant quadrupole moment, and 



1/3 



(15) 



is the orbital velocity. Ryan's formula is valid for any 
spacetime satisfying the assumptions listed above. The 
quadrupole moment of a selected spacetime is therefore 
determined by computing AE explicitly for this space- 
time, and comparing with Eq. (|T^. Because all involved 
quantities are known to be invariant under a coordinate 
transformation, this calculation can be carried out in any 
coordinate system. 

We now go through the steps of calculating AE for the 
specific spacetime of Eqs. (0)-(0). 

Starting from the r component of the geodesic equa- 
tion, using the fact that r and 9 are both constant on 
circular, equatorial orbits, we easily derive the following 
expression for the angular velocity: 



n 



dt 



(16) 



where a comma indicates partial differentiation, and all 
metric functions are evaluated at 6 — tt/2. Reading off 
the metric components from Eq. (Q), substituting the 
asymptotic relations (P)-(p^), taking a cubic root, and 
finally, expanding in powers of a; = (Af/r)^/^, yields 

v = x-^x^~^Xx^ + ^^'+Oix^). (17) 
This series can be inverted to give 

x^v + lv'' + l-xv^ + ^^v^ + 0{v''). (18) 
2 3 4 ^ ' ^ ' 

An expression for E follows from the fact that in 
the coordinates of Eq. (^, £^ = —gtaU°', where = 
(m*, 0, 0, u"^) = u*(l, 0, 0, il) is the particle's four-velocity. 
Using also the normalization condition gapu^u^ = —1, 
we arrive at 



E : 
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(19) 



Going through the same steps as before, we obtain 

E=l-\x^ + \x^-Xx^ + ^^x^ + 0{x''l (20) 



or, after substituting Eq. (|1 
E = l--v^ + -v^~-XV^ 



27-8(7 6 



16 



+ O(w^). 



(21) 



Finally, Eqs. (^ and (|T|) give AE = ~{v/3)dE/dv, 
and substituting Eq. (pi]), we arrive at 



AE^-.^ 



20 



XV 



27 



(22) 



Comparing with the general relation (h4|) confirms that 
J = xM^, as was first stated in Eq. (|11|), and establishes 
that the quadrupole moment is given by Q = qM^, in 
accordance with our original guess. 

We conclude that the dimensionless parameter q ap- 
pearing in Eq. ^ is related to the coordinate invariant 
multipole moment Q by the relation 



qM\ 



(23) 

which is formally identical to Eq. (H). We see that the 
quadrupole moment is indeed determined by isolating the 
"P2(cos0)/r^" part of the relativistic analogue of the 
Newtonian potential. 



III. COMPUTATION OF THE QUADRUPOLE 
MOMENT 

The computation of the quadrupole moment of a neu- 
tron star requires a computer code capable of solving the 
hydrostatic and Einstein field equations for uniformly ro- 
tating mass distributions, under the assumptions of sta- 
tionarity, axial symmetry about the rotation axis, and 
reflection symmetry about the equatorial plane. Along 
with the density profile e(r, 9) and other matter vari- 
ables, the code must compute the metric functions 1^, 
B, to, and a. Such a code was recently written and made 
publicly available by Stergioulas (1995). [See Stergioulas 
& Friedman (1995) for a comparison with other codes.] 
The code, named rns, uses tabulated equations of state 
for neutron-star matter, and is based upon the numeri- 
cal methods of Komatsu, Eriguchi, and Hachisu (1989). 
As a technical point, we may remark that rns uses the 
metric variables p = 2v — hiB and 7 = ln_B instead of v 
and B. 

Once i/(r, 9) has been obtained for a given equation of 
state and for selected values of stellar parameters, such 
as gravitational mass M and angular momentum J, the 
computation of q is in principle straightforward. Indeed, 
Eq. (§) immediately implies 
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-- lim [r/Mf 

2 r— *oo 



i/(r, 61) P2 (cos 6')d cos 6*, (24) 



TABLE I. Properties of equations of state. The entries for 
Mmax are taken from Cook, Shapiro, and Teukolsky (1994). 



and the integral to the right-hand side can easily be 
evaluated numerically. The difficulty with implement- 
ing this method is that when r is much larger than M , 
the "P2(cos6')/r'^" part of v is extremely small compared 
with its spherically symmetric part. The method there- 
fore lacks numerical accuracy. Nevertheless, we have 
found Eq. (p^ ) useful as a check on the results of another 
method of computing g, which we now describe. 

This method was first devised by Salgado et 
al. (1994a), and then independently by Ryan (1997), who 
shows that the quadrupole moment can be straightfor- 
wardly computed as an integral over Sp{r',9'), a com- 
plicated function involving matter variables and metric 
components. This function is defined in Eq. (10) of Ko- 
matsu, Eriguchi, and Hachisu (1989). Ryan's expression 
is 

Q = ^ J Spir\e')r'^P^icos9')d^x'. (25) 

This is the gen eral relativistic analogue of Eq. (|5|). Al- 
though Eq. (|25| ) looks quite simple, many complications 
associated with the nonlinearities of general relativity are 
hidden in the function Sp which, unlike the mass density 
£, has support both inside and outside the star. It is only 
in the Newtonian limit that Sp — 8tt£. Nevertheless, this 
method is easy to implement because rns computes Sp 
and performs similar integrations in order to solve the 
Einstein field equations. It was therefore straightforward 
to make the necessary modifications to the code, and 
compute q = Q/M^ for selected configurations of rotat- 
ing neutron stars. 



IV. RESULTS AND DISCUSSION 

With rns we have constructed numerical models of ro- 
tating neutron stars for four different tabulated equations 
of state, those labeled G, EPS, C, and L. [See Cook, 
Shapiro, and Teukolsky (1994) for a detailed description 
of these equations of state.] For five selected values of 
the gravitational mass (in the interval between 1.0 Mq 
and 1.8 Mq) and for each equation of state, we have con- 
structed models of varying angular momentum J, from 
J = to the Keplerian limit J = Jmax, at which the star's 
angular velocity exceeds that of a test particle in circular, 
equatorial motion just outside the stellar surface. [Of the 
five selected values for the mass, M — lA Mq is probably 
the most relevant, because it is the most commonly ob- 
served for neutron stars in binary systems (Finn 1994).] 
For each of these neutron-star models we have computed 
the quadrupole moment Q using the two methods de- 
scribed in the previous section. 

Table 1 displays some properties of the selected equa- 
tions of state. The first column lists the equations of 



EOS M^ax/Me for J = Xmax 


for M = 1.4Mq 


G 1.36 


0.616 


EPS 1.80 


0.669 


C 1.86 


0.654 


L 2.70 


0.731 


TABLE 11. Quadrupole moment: M 


= 1.0 Mq 


EOS G EOS EPS EOS C 


EOS L 


X <1 X <1 X Q 


X <? 



0.173 


-0.158 


0.182 


-0.289 


0.166 


-0.279 


0.172 


-0.396 


0.232 


-0.275 


0.245 


-0.514 


0.257 


-0.637 


0.269 


-0.960 


0.279 


-0.391 


0.293 


-0.709 


0.325 


-1.028 


0.340 


-1.513 


0.320 


-0.515 


0.336 


-0.921 


0.382 


-1.404 


0.400 


-2.053 


0.358 


-0.642 


0.373 


-1.125 


0.432 


-1.786 


0.452 


-2.578 


0.393 


-0.774 


0.408 


-1.329 


0.477 


-2.164 


0.500 


-3.097 


0.425 


-0.909 


0.440 


-1.540 


0.519 


-2.541 


0.543 


-3.608 


0.460 


-1.067 


0.469 


-1.744 


0.559 


-2.928 


0.584 


-4.118 


0.495 


-1.239 


0.497 


-1.955 


0.596 


-3.319 


0.622 


-4.625 


0.531 


-1.434 


0.525 


-2.167 


0.631 


-3.713 


0.658 


-5.126 


0.566 


-1.635 


0.551 


-2.375 


0.632 


-3.720 


0.693 


-5.629 


0.599 


-1.844 


0.644 


-3.217 











state, in order of increasing stiffness. The second column 
lists Mijiax (in solar masses) for J = 0, the maximum 
value of the gravitational mass for a nonrotating config- 
uration. The third column lists Xmax = cJmax/GM^ for 
M = 1.4 M0. The first row of Table 1 indicates that 
EOS G is extremely soft, with a maximum mass (in the 
absence of rotation) lower than 1.4 Mq. The models con- 
structed for this equation of state therefore belong to the 
supramassive class (Cook, Shapiro, & Teukolsky 1994), 
and are probably not realistic. Furthermore, even rota- 
tion cannot support such a star if its mass is larger than 
1.56 solar masses (Cook, Shapiro, & Teukolsky 1994). 

Our results are presented in Tables 2-6. For each of the 
equations of state, the first column lists x = cJ/GM'^, 
the dimensionless angular-momentum parameter, and 
the second column lists the corresponding value oi q = 
c^Q/G^M^, the dimensionless quadrupole-momcnt pa- 
rameter. The last entry corresponds to a configuration 
of maximum rotation. We recall that neutron stars of 1.6 
and 1.8 solar masses cannot be constructed with EOS G. 
As the number of digits quoted in the tables indicates, 
we believe our computations to be accurate at least to 
one part in 10'^. 

The fact that the quadrupole moments are all negative 
reflects the oblateness of the mass distribution; the rota- 
tion produces a bulge at the equator. For given values 
of M and Xi we find that \q\ increases with the stiffness 
of the equation of state. Thus, \q{G)\ < |g(FPS)| < 
\q{G)\ < |9(L)|. This was to be expected, because a 
stiffer equation of state produces a larger star, and the 
quadrupole moment scales with the square of the star's 
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0.1 0.2 0.3 0.4 0.5 0.6 



X 

FIG. 1. Plots of g as a function of x iov four equations of 
state, for neutron stars of mass M — 1.0 Mq. The data of 
Table 2 is shown with circles, and the solid curves represent 
the fits of Eq. (H). 

TABLE III. Quadrupole moment: M = 1.2 Mq 



EOS G EOS EPS EOS C EOS L 



X 


q 


X 


q 




X 




<3 




X 




q 


0.203 


-0.124 


0.191 


-0.209 





175 


-0 


213 





146 


-0 


246 


0.283 


-0.241 


0.276 


-0.448 





222 


-0 


351 





258 


-0 


683 


0.345 


-0.363 


0.341 


-0.676 





259 


-0 


475 





334 


-1 


120 


0.398 


-0.491 


0.398 


-0.920 





292 


-0 


596 





396 


-1 


546 


0.445 


-0.625 


0.448 


-1.156 





322 


-0 


719 





450 


-1 


975 


0.488 


-0.765 


0.494 


-1.394 





349 


-0 


845 





499 


-2 


395 


0.528 


-0.910 


0.536 


-1.639 





376 


-0 


969 





544 


-2 


809 


0.563 


-1.054 


0.577 


-1.892 





399 


-1 


088 





586 


-3 


218 


0.596 


-1.206 


0.615 


-2.148 





422 


-1 


216 





625 


-3 


623 






0.651 


-2.412 





444 


-1 


343 





663 


-4 


032 













515 


-1 


805 





698 


-4 


439 













579 


-2 


275 





















638 


-2 


768 











TABLE IV. Quadrupole moment: M = 1.4 Mq 



EOS G EOS EPS EOS C EOS L 









n 

y 




A 




n 

y 




A 




y 




A 




y 





337 


-0 


173 





119 


-0 


074 





200 


-0 


199 





124 


-0 


134 





354 


-0 


196 





224 


-0 


216 





269 


-0 


356 





234 


-0 


454 





372 


-0 


223 





292 


-0 


355 





325 


-0 


523 





308 


-0 


764 





393 


-0 


256 





347 


-0 


501 





373 


-0 


690 





368 


-1 


070 





416 


-0 


295 





396 


-0 


655 





417 


-0 


861 





420 


-1 


373 





440 


-0 


341 





440 


-0 


808 





457 


-1 


036 





467 


-1 


672 





466 


-0 


394 





481 


-0 


970 





494 


-1 


211 





510 


-1 


975 





493 


-0 


456 





519 


-1 


134 





529 


-1 


392 





551 


-2 


279 





521 


-0 


526 





556 


-1 


302 





562 


-1 


578 





589 


-2 


587 





549 


-0 


605 





591 


-1 


476 





593 


-1 


768 





626 


-2 


894 





578 


-0 


696 





625 


-1 


659 





623 


-1 


962 





661 


-3 


199 





606 


-0 


795 





658 


-1 


847 





653 


-2 


166 





694 


-3 


507 
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FIG. 3. Same as Fig. 1 but with M = 1.4 Mq. The circles 
represent the data of Table 4. 
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X 

FIG. 2. Same as Fig. 1 but with M = 1.2 Mq. The circles 
represent the data of Table 3. 



TABLE V. Quadrupole moment: M = 1.6 Mq 



EOS FPS EOS C EOS L 



X 


q 


X 


q 


X 


<? 


0.111 


-0.037 


0.130 


-0.066 


0.131 


-0.120 


0.225 


-0.158 


0.223 


-0.184 


0.246 


-0.388 


0.299 


-0.267 


0.286 


-0.289 


0.324 


-0.672 


0.361 


-0.390 


0.340 


-0.409 


0.387 


-0.944 


0.415 


-0.517 


0.386 


-0.528 


0.442 


-1.216 


0.464 


-0.652 


0.429 


-0.658 


0.491 


-1.488 


0.509 


-0.797 


0.468 


-0.793 


0.537 


-1.760 


0.552 


-0.953 


0.506 


-0.937 


0.580 


-2.033 


0.594 


-1.120 


0.543 


-1.089 


0.620 


-2.310 


0.635 


-1.297 


0.578 


-1.248 


0.658 


-2.589 






0.612 


-1.417 


0.695 


-2.872 






0.645 


-1.596 


0.731 


-3.162 
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TABLE VII. Fit parameter 



-1 



-2 



0.1 0.2 0.3 0.^. 



0.5 0.6 0.7 

X 



FIG. 4. Same as Fig. 1 but with PI = 1.6 Mq. The circles 
represent the data of Table 5. 



TABLE VI. Quadrupole moment; M = 1.8 Mq 



EOS EPS 



EOS C 



EOS L 



X 


q 


X 


9 


X 


q 


0.113 


-0.024 


0.112 


-0.030 


0.136 


-0.109 


0.169 


-0.052 


0.191 


-0.081 


0.250 


-0.326 


0.222 


-0.091 


0.253 


-0.152 


0.326 


-0.549 


0.274 


-0.143 


0.304 


-0.222 


0.388 


-0.768 


0.323 


-0.203 


0.353 


-0.306 


0.443 


-0.990 


0.372 


-0.278 


0.399 


-0.396 


0.491 


-1.214 


0.420 


-0.365 


0.442 


-0.496 
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FIG. 5. Same as Fig. 1 but with M = 1.8 A/q. The circles 
represent the data of Table 6. 
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For fixed mass and equation 
on X is well reproduced with 


of state, the depen- 
a simple quadratic 




q - 


i -a{M, EOS) 




(26) 


or 




a(M, EOS) J2 
c2 M 




(27) 



with the parameter a depending on the mass and the 
equation of state. The best-fit values are presented in 
Table 7, which shows that for a fixed equation of state, a 
decreases with increasing mass. (The relation is not well 
reproduced by simple fitting formulas, such as linear and 
power-law relationships.) The accuracy of these fits is 
demonstrated in Figs. 1-5, which display g as a function 
of X for each of the selected equations of state. The 
figures show the data points of Table 2-6 as well as the 
fits. We see that Eq. (|2^) reproduces the data quite well. 

The high degree of accuracy of the formula q ~ —ax^ 
is intriguing, because it holds not only for slow rotations, 
where such a quadratic relation is to be expected, but 
also for fast rotations, where such a relation is expected 
to fail. Here, our expectation is based on the theory of 
uniformly-rotating, constant-density fluid configurations 
in Newtonian gravity (Chandrasekhar 1969). For such 
configurations, called Maclaurin spheriods, the relation 
between q and x is known exactly. For slow rotation, the 
relation is quadratic: q ~ ~{25/8){c^Sp/GM) x^ , where 
Sp is the polar radius. For fast rotations, however, the 
relation deviates strongly from a quadratic. For typical 
neutron-star parameters (M = 1.4 Mq and Sp — 15 km), 
the deviation becomes significant when x ^ 0.4, which is 
smaller than the maximum value (Ymax > 0.6). There- 
fore, the reliability of the relation (|26| ) in the complete 
interval < x < Xmax is quite remarkable. In this regard, 
we may note that the slow-rotation approximation was 
shown by Weber and Glendenning (1992) to give rather 
accurate results in numerical models of rotating neutron 
stars. We also note that the relation q = —y^ is known 
to hold exactly for rotating black holes (Thorne 1980). 
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